Local quantum control of Heisenberg spin chains 
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Motivated by some recent results of quantum control tlieory, we discuss the feasibility of local 
operator control in arrays of interacting qubits modeled as isotropic Heisenberg spin chains. Acting 
on one of the end spins, we aim at finding piecewise-constant control pulses that lead to optimal 
fidelities for a chosen set of quantum gates. We analyze the robustness of the obtained results for 
the gate fidelities to random errors in the control fields, finding that with faster switching between 
piecewise-constant controls the system is less susceptible to these errors. The observed behavior falls 
into a generic class of physical phenomena that are related to a competition between resonance- and 
relaxation-type behavior, exemplified by motional narrowing in NMR experiments. Finally, we 
discuss how the obtained optimal gate fidelities are altered when the corresponding rapidly-varying 
piecewise-constant control fields are smoothened through spectral filtering. 



PACS numbers: 03.67.Hk, 03.67.Lx, 75.10.Pq 



I. INTRODUCTION 

The ability to manipulate the state and/or the dynam- 
ics of complex quantum systems is one of the main ob- 
jectives of quantum information science [1| . In this con- 
text, quantum control i2| has been successfully utilized 
to address two general classes of issues. In state- selective 
control, the main question is how to steer a physical sys- 
tem from an initial reference state to a desired final state. 
In operator (state-independent) control, one seeks to im- 
plement a predetermined unitary transformation (target 
quantum gate) irrespective of the initial state of the sys- 
tem, which is often unknown. The rigorous mathematical 
foundations of the field, based on the notion of controlla- 
bility and framed using Lie-algebraic concepts, have long 
been known [3[. Recent quantum-control studies have 
been focused around two sets of questions. 

Firstly, even when a system is fully controllable in prin- 
ciple, it is of interest to know how to control it most effi- 
ciently, taking into account various practical constraints. 
Issues of this type are not likely to yield universal an- 
swers, except for well-understood general topological fea- 
tures of optimal-control landscapes jj] . Secondly, it is de- 
sirable to know whether the system can be partly or fully 
controlled by acting only on a small subsystem. This is 
the main idea behind the local- control (minimal actua- 
tion) approach. Such an approach is applicable only in 
interacting systems like coupled spin chains, where the 
interaction will effectively turn the local control applied 
to the small subsystem to a global one. Uses of such 
"always-on" interactions in quantum information pro- 
cessinginclude, e.g., systems which can serve as data 
buses enabling state- p-Q and entanglement trans- 
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fer [T3|. Several notable results have recently been re- 
ported pT| - [l4| . For instance, controlling only one of 
the end spins of an XXZ-Heisenberg spin chain ensures 
full controllability of the chain [I4I . In addition, univer- 
sal quantum computation with the XY spin chain can 
be effected by controlling two spins at the end of the 
chain [l3[ . Finally, a magnetic field in the z-direction act- 
ing only on a single spin suffices for generating perfect en- 
tanglement between the ends of a Heisenberg chain |14[ . 

Apart from being interesting from a conceptual point 
of view, the local-control approach also has some prac- 
tical importance. In addition to its easier implementa- 
tion, lowering the number of control actuators reduces 
the effects of decoherence. This is of crucial importance 
in many quantum-computer architectures, such as super- 
conducting qubits [1^ [ig . 

In thispaper, motivated in part by the rigorous results 
of Ref. 13, we consider an implementation of local opera- 
tor control in arrays of interacting qubits modeled as spin 
chains with isotropic Heisenberg interaction. Acting only 
on one of the end spins, we determine piecewise-constant 
control fields resulting in the highest possible fidelities 
for a chosen set of quantum gates. By treating the am- 
plitudes of these control fields as optimization variables, 
we find the optimal fidelities for chains with up to four 
spins 112] . While our work focuses on Heisenberg-coupled 
spin chains, its results are relevant for any physical im- 
plementation of interacting qubit arrays. 

While some elements of the quantum control of spin 
chains have already been studied [l^ - even in the open- 
system scenario [l9| - we discuss in detail some aspects 
that have so far not received due attention. In partic- 
ular, we carry out a sensitivity analysis, i.e., consider 
the robustness of the obtained results for the gate fideli- 
ties with respect to random errors in the control fields. 
Our analysis shows that with faster switching between 
piecewise-constant controls, the system is less suscepti- 



ble to these errors. Importantly, we explain this behavior 
by making a link to a class of phenomena exemplified by 
motional narrowing in NMR experiments [20| • To bridge 
the gap between theoretical considerations and future ex- 
perimental implementations of the system under study, 
we discuss how the optimal fidelities are affected when 
the control pulses are smoothened by eliminating high- 
frequency components in their Fourier spectra (spectral 
filtering). We establish some qualitative criteria regard- 
ing the performance of such an approach. 

The paper is organized as follows. To set the stage, in 
Sec. |TT] we introduce the system and our control objec- 
tives. In Sec.|nT]we present detailed results for the opti- 
mal gate fidelities and the corresponding control fields in 
chains with three and four spins. We also discuss mini- 
mal action times for the control fields needed to realize 
certain quantum gates. Section IIVI is set aside for the 
sensitivity analysis, while Sec. FVl deals with the effects of 
spectral filtering of the optimal piecewise-constant con- 
trol fields on the resulting gate fidelities. We conclude 
with a brief summary of the paper and some general re- 
marks in Sec. IVII Some numerical details are described 
in the Appendix. 



II. SYSTEM AND METHOD 

A. Hamiltonian and control objectives 

We consider an isotropic Heisenberg spin- 1/2 chain of 
length Ns , with control fields acting on the first spin only. 
It is governed by the total Hamiltonian 

H{t)=Ho + H,it), (!) 

where Hq is the Heisenberg Hamiltonian 

Ho = J y ^ {SixSi+l^x + SiySiJ^l^y + SizSi+l^z) 7 (2) 



while 



Hc{t) = hx{t)Six + hy{t)Sly 



(3) 



is the Zeeman-like control Hamiltonian. The time depen- 
dence of the control fields hx{t) and hy{t) will be speci- 
fied shortly. For convenience, throughout this paper we 
set S = 1. As a consequence, all frequencies and control- 
field amplitudes can be expressed in units of the coupling 
strength J, and all times in units of 1/J. 

Regarding the choice of Hamiltonian in Eq. ([2]) two 
remarks are in order. Firstly, while the algebraic re- 
sult of Ref. 1I2! is even applicable to the more general 
(anisotropic) AAZ-Heisenberg Hamiltonian, for simplic- 
ity we here discuss the isotropic case only. Secondly, 
whether the spin chain is ferromagnetic {J < 0) or anti- 
ferromagnetic ( J > 0) is unessential for our present pur- 
poses, because we are concerned with operator control. 
For concreteness we assume that J > 0. 



In Ref. |l2|, using the graph infection criterion, it was 
shown that controlling the x and y components of the 
field acting on the first spin guarantees the complete con- 
trollability of a Heisenberg chain. Moreover, there are 
unitary transformations that require even smaller degree 
of control, i.e., a control field only in one direction. For 
instance, to achieve a spin-flip operation X on the last 
spin of the chain, where X is the Pauli matrix, one needs 
only a control field in the x direction. The corresponding 
dynamical Lie algebra Cx, generated by {—iHo, —iSix}, 
is a subalgebra of su{d). The action of the A^^-qubit gate 
Ajv^ , which flips the last spin in the Heisenberg chain, is 
defined by 



X 



Ns 



l(g)l« 



1® A 



(4) 



To prove the reachability of X^^, i.e., that Xn^ e e^"" 
[where e^" is the connected Lie subgroup of SU{d) with 
Lie algebra Cx], we have to show that there exists an 
element A of the dynamical Lie algebra Cx such that 
Ajv^ = e^. By calculating the (repeated) commutators 
of the operators which generate the algebra, we find that 
AjVj, is an element of Cx ■ Since A^r^ is both unitary and 
Hermitian (hence A^ = 1), the operator A — —i^XN_^, 
an element of Cx, has the property that e = — iA^r^. 
This concludes the proof that the x field is sufficient to 
reach Xn^ . More generally, any unitary U, which is also 
Hermitian and for which —iU belongs to the dynamical 
Lie algebra C, is an element of the reachable set e^. 

In the following, our control objective is the realization 
of concrete quantum gates. We discuss gates that require 
control fields in the x and y directions, as well as those 
that entail only an x control field. Apart from spin-flip 
(NOT) gates, we implement entangling two-qubit gates 
such as 



CNOT 



JVs 



1® 1 



1 (E) CNOT 



(5) 



which performs the controUed-NOT (CNOT) operation 
on the last two qubits in the chain. Another example 
is square root of SWAP (vSWAP). In operator control, 
the figure of merit in the realization of these gates is the 
gate fidelity 



F(i) = ^|tr[C/t(i)C/t,,get]| 



(6) 



where U{t) is the time-evolution operator of the system 
at time t and C/targct stands for the desired quantum gate. 



B. Time-evolution for piecew^ise-constant control 
fields 

For many classes of systems, the implementation of 
complicated time-dependent potentials is rather difficult. 
In what follows, we resort to simple piecewise-constant 
controls. Importantly, we retain the full Hilbert space of 
the system, unlike some previous studies that make use of 



the single-excitation subspace. [22] This puts constraints 
on the system size. 

Assume that we want to achieve an arbitrary target 
unitary at a time t — tf. At i = we apply an x control 
pulse to the first spin of the chain with amplitude h^i 
which is constant throughout the pulse duration T. That 
is, the system evolves under the action of the Hamilto- 
nian H^^i = Hq + hx.iSix- Then we apply a y control 
pulse with amplitude hy^i during the second time inter- 
val of length T, whereby the system is governed by the 
Hamiltonian Hy i = Hq + hy iSiy. This sequence of al- 
ternate X and y control pulses is repeated until Nt pulses 
have been completed at the time tf = NfT. The full 
time-evolution operator can be expressed as 



Uitf) — Uy^N^/2 ■ Ux,Nt/2 ■ ■■■ ■ Uy.l ■ Ux,l 



(7) 



where Ux,i = e"*^"''-'" and Uy^i = e^*^"''^ are the re- 
spective time-evolution operators corresponding to H^^i 
and Hy^i. We evaluate Ux.i and Uy^i using their spectral 
forms. If the desired gate is achievable by control of the 
X field only, it is sufficient to apply in each time interval 
an X field of variable amplitude. 



III. THREE- AND FOUR-SPIN CHAINS 

For the numerical maximization of the fidelity [Eq. ([5])] 
with respect to the Nt field amplitudes, we use a quasi- 
Newton method developed by Broyden, Fletcher, Gold- 
farb, and Shanno (BFGS-algorithm) [24]. We first choose 
an initial guess for the control-field amplitudes. The al- 
gorithm then generates iteratively new sequences of field 
amplitudes such that at each iteration point the fidelity 
is increased, and terminates as soon as the desired accu- 
racy is reached. This procedure ensures the convergence 
to a local maximum, but of course does not guarantee 
the convergence to a globally-optimal sequence of con- 
trol field amplitudes. 

We perform maximizations with varying number Nt 
and durations T of pulses, and hence different total evo- 
lution times tf = NtT. It is important to stress that, 
instead of fixing the pulse durations and maximizing over 
the field amplitudes, we could as well keep the amplitudes 
constant and treat the pulse durations as free control pa- 
rameters [111 ? ]. Yet, we choose to maximize over the 
field amplitudes as this approach allows us to fix easily 
the total evolution time and estimate its minimal value 
needed to implement the desired gate. Although control- 
lability implies the existence of a control that can enable 
realization of a desired quantum gate, it does not provide 
any information about the minimal time over which this 
realization is possible. 

We now give explicit examples of control sequences 
that implement all of our chosen gates in chains of three 
and four spins, among them entangling gates which are 
essential for universal quantum computation. Two opti- 
mal sequences of control pulses implementing the spin- 
flip gate A'tVj for a three-spin chain (Ng = 3) with x 
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FIG. 1: (Color online) Examples of optimal control fields for 
the X3 gate, corresponding to: (a) Aft = 70 and T = 0.429, 
with the resulting optimal fidelity of 1 — 10~*, and (b) Nt = 
70 and T — 0.243 (minimal total time tf — 17), with the 
resulting optimal fidelity of 1 — 10~ . 



and y pulses are shown in Fig. [T] [the optimal held in 
Fig. [Ijb) corresponds to the minimal time needed]. For 
given total time longer than the minimal time, we can 
generate piecewise-constant control fields with fidelities 
arbitrarily close to unity by increasing the switching rate. 
For instance, in the case of implementing the gate X3 by 
control of the x and y fields, for a total time of tf = 30 
and Nt = 10, 20, 30, 40, 50, 60, 70, we respectively obtain 
F = 0.536, 0.764, 0.904, 0.964, 0.992, 0.999, 1 - 10"^ 

We give here an estimate of the minimal time for the 
realization of the gate X^^ (operating on three- and four- 
spin chains) by controlling the x and y field. The im- 
plementation of this unitary is possible for A^s = 3 and 
Ns = i within the respective times tf = 17 and tf = 70, 
which are lower bounds on the evolution time needed. 
Unlike in Ref. [ij, the size of the dynamical Lie algre- 
bra scales exponentially with the number of qubits. We 
therefore expect that the minimal time grows rapidly 
with the chain length, and this is a limiting factor for 
extending this control procedure. 

For each fixed chain length, the minimal times for re- 
alizing the spin-flip and CNOT gates are similar. This 
result seems somewhat related to the main conclusion of 
Ref. [2^ where the CNOT gate is implemented for two 
qubits coupled through a variety of interactions {XY, 
Heisenberg, Ising) and the total gate time does not de- 
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FIG. 2; (Color online) Average fidelity versus half-width (S) 
for random noise affected optimal fields corresponding to (a) 
Xa gate, and (b) CNOT3 gate, both implemented by alternate 
X and y controls with a total evolution time t/ = 30. 



pend significantly on the choice of interaction as long as 
it is not of pure Ising type. 

Recalling that gates like X^^ require only an x control 
field, it is interesting to compare the minimal times for 
their implementation depending on the degree of control. 
We find that with control of the x field only, the respec- 
tive gates can be achieved in approximately the same 
time as with control of both the x and y fields. 



IV. SENSITIVITY ANALYSIS 
A. Preliminary considerations 

In the following, we analyze the sensitivity of the fi- 
delity to random errors in the control fields. To this end, 
we add random numbers from a uniform distribution of 
half-width 5 to the optimal control-field amplitudes. For 
given S, we generate a few hundred sequences of ran- 
dom numbers. In this way we obtain a large sample of 
N ~ 1000 control fields affected by random noise, for 
which we recalculate the fidelity. 

In the following, we discuss the behavior of the average 
fidelity F = J2i=i ^i/-^j where the Fi are fidelities cal- 
culated for specific realizations of the random field. We 
note that the standard deviation <Tp increases with S: for 
6 = 0.01 it is found to be of the order of 10^^, while 
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FIG. 3; (Color online) Average fidelity versus half-width (S) 
for random noise affected optimal fields corresponding to the 
implementation of the X3 gate by control of the x field only: 
(a) fixed total evolution time (i/ = 25) and (b) fixed number 
of pulses {Nt — 70). The insets show the same curves for 
large values of the parameter S. 



for S = 1.0 we obtain ap ^ 0.1. Figures [5] and [3] illus- 
trate how the strength of the randomness in the optimal 
field amplitudes affects the average fidelity. While Fig. [2] 
shows the results for the gates X3 and CNOT3 both re- 
alized by a sequence of alternate x and y pulses, Fig. [3] 
refers to the implementation of X3 by control of the x 
field only. We see that F is less susceptible to random 
noise for shorter pulse durations T. A physical under- 
standing of this phenomenon is provided in Sec. I IV CI 

Starting from optimal fields, the shape of the fidelity 
decay with increasing S is dependent on the number of 
control pulses Nt and their length T, but not on the 
pulse amplitudes. The saturation regime of the average 
fidelity, indicating full randomization, sets in for S > J. 
Importantly, the saturation of the fidelity is an intrinsic 
property of our system. If the system is completely con- 
trollable, the saturation value is universal, i.e., indepen- 
dent of the concrete shape of the control field, provided 
the evolution time tf and the number of pulses Nt al- 
low the generation of any unitary (illustrated in Fig. [2] 
for tf — 30). In contrast, if the dynamical Lie algebra 
of the system is a proper subalgebra of su{d), for differ- 
ent evolution times tf the average fidelity may saturate 
at different values, and the saturation may require much 
stronger random fields in comparison to the completely 
controllable system (for illustration, see Fig. |3]). 



B. Fidelity saturation 

To understand the value of the fidehty at saturation, 
we invoke the notion of the average gate fidehty [2g . If 
e is a trace-preserving quantum operation pertaining to 
a quantum system with d-dimensional Hilbert space and 
U a certain quantum gate, the average gate fidehty 



Fie,U) 



E-litr[c/C/jC/t£(C/, 



d2(d+l) 



(8) 



where {Uj} is an orthonormal unitary-operator basis of 
the complex vector space C''^'', quantifies how well e 
approximates U. If e implements the gate U perfectly, 
this mapping is given hy e : p —> e{p) = UpW , implying 
that F{e, U) = I; otherwise, the implementation is noisy 
and, consequently, F{e,U) < 1. An equivalent form of 
Eq. ([5]) - which involves the traceless generators |r,j of 
SU{d) [recah that tr{TjTk) = <5jfe/2, see e.g. Ref. 
reads [27j: 



21| 



,U) 



cf-l 



d{d + l) ^ 



^ tr [UT,U^e{T,)\ 



(9) 



Equations ^ and (O are both derived by integrating 
over the uniform Haar measure [26|, hence requiring 
a uniform starting distribution of states and/or opera- 
tors. Since the quantum control system governed by the 
Hamiltonian H{t) = Hq -\- hx{t)Six + hy{t)Siy allows to 
generate any element of SU{d), Eq. ([9]), which is derived 
explicitly in terms of the SU{d) group generators, pro- 
vides a suitable expression for calculating the average 
gate fidelity. Assuming that the condition of full random- 
ization is fulfilled, the action of the quantum operation e 
is given by 



e : p-^ e{p) 



1 
1 



(10) 



for every p, that is, an arbitrary (pure) state is mapped 
onto a maximally-mixed state. We recall that an arbi- 
trary density matrix can be written as 



1 



= 1 + ^2"^^^ («^-^^) 



(11) 



Since s{p) = t/d for an arbitrary p, by assumption, it 
follows by linearity and trace-preserving that e{Tj) = 0, 
for every j = 1, . . . , d^ — 1. Hence, the sum in Eq. © 
evaluates to 0, implying that the average gate fidelity for 
the mapping in Eq. (|10p is given by 



Fie,U) 



1 



fl2) 



Note that Eq. p^ is applicable only if the randomness 
in the control field amplitudes allows the uniform gener- 
ation of any unitary contained in SU{d). For example. 



in the three-spin (d = 8) case we expect the average 
gate fidelity to saturate at a value of 0.125. This seems 
to be numerically corroborated by Fig. [5] and indicates 
that the system undergoes full randomization. Also in 
a four-spin chain {d — 16), provided that full random- 
ization occurs, the found saturation value fits well with 
l/d = 0.0625 (not illustrated here). It should be em- 
phasized that this regime of full randomization is of no 
relevance for practical realizations of control, however, it 
provides an important consistency check of our numerical 
results. 

We now address the observation that there exists no 
universal saturation value in the case of control by x 
pulses only, see Fig. [S] This type of control does not 
allow universal quantum computation, i.e., the reachable 
set of unitary operators is reduced to a Lie subgroup of 
SU{d). Since Eq. ([5]) is only valid if the reachable set 
equals the full SU{d), the reasoning leading to Eq. p2)) 
does not apply to such a subgroup. The accessible part of 
the reachable set depends on the evolution time tf . The 
control fields of Fig. [3jb) differ in tf which may result in 
different accessible parts of the reachable sets and, con- 
sequently, different saturation values. The fact that in 
Fig. El^a), where all fields have the same tf, the average 
gate fidelity saturates to the same value corroborates this 
explanation. 



C. Physical interpretation 

As illustrated in Figs. ^ and |31 the average fidelity is 
less susceptible to random errors for shorter time inter- 
vals T. This is a special case of a more general physical 
situation, namely a competition between resonance- and 
relaxation- type behavior [28]. A typical example is the 
phenomenon of motional narrowing of the linewidth in 
NMR experiments [20J . 

This generic class of phenomena is discussed based 
on models for a randomly-interrupted deterministic mo- 
tion [2^: a system is subjected to a random event ( e.g., 
switching between magnetic fields Bq and —Bq in NMR) 
with the switching rate A, while in between two random 
events it evolves deterministically. On the whole, the sys- 
tem is then evolving piecewise-deterministically and its 
qualitative behavior depends on the relative magnitude 
of the switching rate A and the average (absolute) exter- 
nal field amplitude {Bq in the example above). Note that 
with units chosen such that h — I both quantities have 
the same physical dimensions. 

In our system, for smaller time intervals T the fidelity 
in the presence of random errors is indeed closer to the 
intrinsic optimal values. Interestingly, A — T^^ ~ J in 
the optimal cases, while the average (absolute) control- 
field amplitudes (that are independent of S because of the 
symmetry of the probability distribution) are of the same 
order of magnitude. This explains the dependence of F 
on S: the fidelity varies significantly with S and saturates 
at a value around 0.125 in the case of a completely con- 



trollable three-spin chain (0.0625 in the four-spin chain). 
In contrast to that, for A much bigger than the average 
control-ampHtudes, the fidehties would remain close to 1 
even for large S. 



V. SPECTRAL FILTERING 

Our choice of piecewise-constant control fields is, at 
least partly, a matter of mathematical convenience. In re- 
alistic implementations of quantum control various con- 
straints may come into play, the most important ones 
pertaining to the complexity of the frequency spectrum 
of permissible control fields. In control experiments that 
make use of an external magnetic field, for instance, such 
constraints are related to bandwidth and slew-rate limi- 
tations of the magnetic coils and drivers [23]. Therefore, 
it is necessary to subject the optimal control fields to 
spectral filtering in order to make contact with experi- 
mental realizations |3(j |. 

In what follows, we aim at finding control fields that 
are smoother than the optimal piecewise-constant ones, 
at the same time retaining high fidelities for our target 
quantum gates. To this end, we put constraints on the 
frequency spectra of the control fields hj(t) {j — x, y) by 
means of frequency filter functions. After operating with 
a filter function /(w) on the Fourier transforms J^[hj] of 
the optimal fields, we switch back to the time domain 
and calculate the filtered control fields hj{t) via inverse 
Fourier transformation 311: 



h,{t) 



^-i[/HJ-[/i,]H](t) ij = x,y). (13) 



The power spectrum of a typical optimal piecewise- 
constant control field is depicted in Fig. ID 

Based on the obtained smoothened control fields, we 
calculate the corresponding fidelities for different quan- 
tum gates, using a product-formula approximation (for 
details, see the Appendix). 

We first consider an ideal low-pass filter which removes 
frequencies outside the interval [^loq^loq]: f{uj) = 9{uj + 
u)q)~9{uj—ujo). Using the general prescription in Eq. (|13p. 
we obtain 



Nt/2 



^ 1 " 

hx{t) = - 2J hx,n[a2n-l{t) - a2n-2{t)\ 



n=l 

Nt/2 



(14) 



^^y{i) = - X! ^y:"['^2n(t) - a2n-l{t)] 



where am{t) = 'A\\i^o{mT ~ t)\ and Si(a::) = J„ (sint/t)dt. 
Two examples of such smoothened pulses for the X^ 
and CNOT3 gates, corresponding to a fidelity of 0.9, are 
shown in Figs. [5] and IH] 

We also consider Gaussian filters with center frequen- 
cies ±uic- /(w) = exp[-7(a; - ujcf] + exp[-7(a; -I- u)c)\ 
where 7 > determines the Gaussian full-width-at-half- 
maximum FWHM = 2\J\a2/^. As the power spectra 




FIG. 4: (Color online) The power spectra corresponding to 
the optimal x [(a)] and y [(b)] control fields shown in Fig.[lla). 



2 



-2 



n ' r 



(a) 



■F=1-10, optimal 
■ F=0.9, (B„/J=7:/2 



id^^ 



LJi \ ^ \ ^ L 




5 10 15 20 25 30 35 
time (units of 1/J) 
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field which realizes the X3 gate by alternate x [(a)] and y [(b)] 
pulses, compared to the low-pass filtered counterpart (cut-ofi^ 
frequency ujq/ J = 7r/2) with a fidelity of 0.9. The optimal 
field corresponds to Nt = 70 and T = 0.5. 
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FIG. 6: (Color online) Optimal piecewise-constant control 
field which realizes the CNOT3 gate by alternate x [(a)] and 
y [(b)] pulses, compared to the low-pass filtered counterpart 
(cut-off frequency luq/J = 7r/2) with a fidelity of 0.9. The 
optimal field corresponds to Nt = 70 and T = 1.0. 

of the optimal control fields (Fig. 0]) show, the dominant 
frequencies are located around a; = 0. Hence we choose 
ujc ~ and vary the Gaussian width. We determine 
the smoothened control fields by making use of the iden- 
tity [34I 



_„2 2sm(aa;) , tt / a 
X 2 \2p 



where erf is the error function. We obtain 



(15) 



Nt/2 



hx{t) = 7^'^ K,n[b2n-l{t) - fo2«-2(i)] 



n=l 

Nt/2 



(16) 



^yW = 9 X! K-yA^2n{t) - b2n-l{t)\ 



where 6™(i) = erf [(mT - t)/(2^)] . The resuhs ob- 
tained through Gaussian filtering resemble the ideal low- 
pass ones. In Fig. [7] we show the optimal sequence of con- 
trols realizing the X^ gate which is now smoothened by 
applying a Gaussian filter, retaining, however, a fidelity 
of 0.9 (compare the corresponding low-pass filtered pulse 
in Fig. El). 

The behavior of the fidelity versus the cut-off (low- 
pass) or the width (Gaussian filtering) turns out to be 
dependent on the strength and pulse shape of the opti- 
mal control fields. Figure |H] illustrates this dependence in 




5 10 15 20 25 30 35 
time (units of 1/J) 

FIG. 7: (Color online) Optimal piecewise-constant control 
field which realizes the X3 gate by alternate x [(a)] and y 
[(b)] pulses, compared to the Gaussian filtered counterpart 
(FWHM = 4.3) with a fidelity of 0.9. The optimal control 
field corresponds to Nt = 70 and T = 0.5 as in the case of 
low-pass filtering shown in Fig. (5] 



the case of ideal low-pass filtering which implement the 
X3 gate. The curves in Fig.ISJa) reveal that the fidelity of 
weaker optimal control fields is less affected by the low- 
pass filtering. The average absolute field amplitude spec- 
ified in the legend has been calculated by averaging the 
absolute values of the pulse amplitudes over the total con- 
trol time tf. In contrast, the optimal fields corresponding 
to the fidelity curves shown in Fig.Hl^b) are characterized 
by a similar average absolute amplitude but differ in the 
variance of the single absolute pulse amplitudes defined 
with respect to the average absolute amplitude. We see 
that the fidelity decays faster if the variance is larger, 
i.e., if the pulse amplitudes of the field are varying more 
rapidly. Analyzing in Fig. [5] the cut-offs at which the 
fidelity decays, we estimate that frequencies of up to at 
least twice the average absolute field amplitude are re- 
quired to retain fidelities of 0.9 or larger. 

To improve the results from the filtering of optimal 
control fields, we tried to iterate the procedure: The fil- 
tered fields are discretized and used as initial guess for 
the optimization algorithm. The optimal fields produced 
by iteration are then filtered again. In so doing we in- 
tended to generate pulses with high fidelity and even 
lower frequency than in the "ordinary filtering" discussed 
before. However, this approach does not yield significant 
improvements. 
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VI. SUMMARY AND CONCLUSIONS 

Recent quantum-control studies have shown that an ar- 
ray of qubits coupled by nearest-neighbor Heisenberg in- 
teractions can be universally operator-controlled by act- 
ing only on the first qubit. Since these results only imply 
the existence of a control sequence but do not provide 
a way to construct them, we have investigated the fea- 
sibility of controlling a spin chain with Heisenberg in- 
teraction by applying a time-dependent magnetic field 
to the first spin in the chain. We have explicitly deter- 
mined piecewise-constant control fields for several quan- 
tum gates for three- and four-spin chains. By increasing 
the number of control pulses within a fixed total time, 
fidelities arbitrarily close to 1 can be achieved. 

We have also studied the sensitivity of the fidelity to 
random errors in the control fields. Our analysis shows 
that the average fidelity is less susceptible to random 
perturbations if the durations of the single control pulses 
are reduced. This behavior is related to a generic class of 
phenomena exemplified by motional narrowing in NMR. 
We have also examined the intrinsic saturation of the av- 
erage fidelity, being universal for complete controllability 
and occurring if the strength of the acting random field is 



large enough. Finally, to make contact with experimen- 
tal realizations, we have used spectral filtering to obtain 
smoothened control pulses retaining fidelities of 0.9. 

Our study can be extended to include more compli- 
cated (e.g., continuously varying) control pulses. Such 
cases would require the use of some advanced methods 
of the optimal control theory, e.g., the Krotov algorithm, 
which was shown to be capable of reaching the quantum 
speed limit Q. 

The statement that universal control of a qubit array 
is possible by acting on the first qubit is true in princi- 
ple for arbitrarily long arrays. However, the complexity 
of finding the control sequence will make such a proce- 
dure impractical for Heisenberg chains. We expect our 
results to be of practical significance for relatively small 
systems of a few qubits, like those realized experimentally 
right now. Having fewer control lines will increase the co- 
herence time of the qubits that are indirectly controlled 
by their interaction with the neighboring qubits. In the 
present paper, we have studied isotropic Heisenberg in- 
teractions, but the procedure also works for XXZ-type 
interactions. We therefore expect that our work will fa- 
cilitate the control of quantum information devices with 
relatively few qubits. 




FIG. 8: (Color online) Fidelity (F) versus cut-ofl frequency 
(a;o) for different low-pass filtered optimal control fields, all 
of them implementing the X3 gate with 70 pulses of dura- 
tions 1.5 alternately applied in the x and y directions: (a) 
The optimal control fields with average amplitudes indicated 
and variances ranging between 1.1 and 1.3. (b) The optimal 
control fields with average amplitudes between 1.9 and 2.1, 
and variances specified in the legend. 
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Appendix A: Fidelity calculation for smoothened 
control fields 

In order to determine the fidelity corresponding to the 
smoothened control fields, we have to calculate the time- 
evolution operator U{t) governed by the Hamiltonian 
H{t) — Hq + hx{t)Six + hy{t)Siy. One possible approach 
to this problem is solving the equation of motion for U{t), 
which can be written as a system of (P coupled first-order 
differential equations 



ij^Ukiit) 






(fc,/ = l,...,d), (Al) 



with Uki (0) — Ski ■ While Runge-Kutta-type methods are 
commonly utilized for solving such systems, they are not 
so convenient in the case at hand where it is essential 
to preserve the unitarity of the time-evolution operator. 
We therefore determine U{tf) using a product- formula 
approximation which manifestly preserves unitarity and 
essentially amounts to a discretization in time. 

We divide each time interval of length T into rrirp 
steps of length t = r/m^,; the total number of steps 
is denoted by mf = Ntmrp. The approximation hinges 



on the assumption that during each time step the total 
Hamihonian of the system remains constant, i.e., time- 
independent. In other words 

H{t) = i/f'^) =Ho + K{kT)Si^ + hy{kT)Siy (A2) 

for kr < t < {k + l)r, with k = 0, . . . , m/ — 1. Hence 
the evolution of the system during the (fc + l)-th interval 



is described by U^''^ (r) 



- .-«-&<" 



Using the semigroup 



property of time-evolution operators, U{tf) can be ex- 
pressed as the product 



where each of the operators e~*^ "^ is found using the 
spectral representation (recall Sec. IIIBj) . This approx- 
imation becomes progressively more accurate with de- 
creasing r, so that the time-evolution operator can be 
computed to the required accuracy. The unitarity of the 
time-evolution operator is preserved by construction for 
an arbitrary number of time steps, hence the method is 
unconditionally stable |33| . 



U{tf) 



-,-iH' 



e-^^*"^•e-^^'"'^ (A3) 
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